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I. ABSTRACT 



III. EXPERIMENTAL SETUP 
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We study liquid-liquid dispersions in a turbulent 
Taylor-Couette flow, produced between two counterro- 
tating coaxial cylinders. In pure Water and in counter- 
rotation, Reynolds numbers up to 1.4 x 10^ are reached. 
The liquids we use are a low-viscous Oil and pure Wa- 
ter or a Sodium Iodide solution with a refractive index 
matched to that of Oil, in order to get transparent dis- 
persions. We first characterize the single-phase flow, in 
terms of threshold for transition to turbulence, scaling of 
the torque and measurements of the mean flow and of the 
Reynolds stress by stereoscopic PIV. We then study the 
increase of the dissipation in the two-phase flows and find 
that the torque per unit mass can be twice the torque for 
a single-phase flow. Long-time behaviours are also re- 
ported. 



II. INTRODUCTION 

Liquid-liquid dispersions are encountered for instance 
in extraction or chemical engineering when contact be- 
tween two liquid phases is needed. Without surfactants, 
when two immiscible fluids are mechanically agitated, 
a dispersed state resulting from a dynamical equilib- 
rium between break-up and coalescence of drops can be 
reached. Their modelling in turbulent flows is still lim- 
ited (Portela and Olicmans, 2006). 

Some open questions we address are first the poten- 
tial increase of turbulent dissipation due to the disper- 
sion, the modification of the turbulent transport inside 
the fiow, and the droplet size distribution, governed by 
the breakup and coalescence rates. Finally, some con- 
centrated liquid-liquid flows are also known to experience 
hysteretic phase inversion (Piela et al., 2006). 

To answer these questions, we have chosen a model 
shear-flow: a Taylor-Couette flow where the fluids are 
stirred in the gap between two counterrotating coax- 
ial cylinders. We use low- viscous fluids and first check 
that the single-phase flows are strongly turbulent. We 
measure the wall shear stress, and use Refractive-Index 
matched fluids in order to measure the phases reparti- 
tion. 



The flow is produced between two coaxial cylinders. 
The inner one is of radius = 110 ± 0.05mm, and the 
outer one of radius = 120 ± 0.05mm, which gives a 
gap ratio r] = ri/vo = 0.917. The length of the cylinders 
is L = 220mm, i.e. the axial aspect ratio is L/{ro — 
Ti) = 22. The cylinders axis is vertical in order to keep 
axisymmetric conditions. The system is closed at both 
ends, with top and bottom lids rotating with the outer 
cylinder. There is a gap of 1.5mm between these lids 
and the flat ends of the inner cylinder. The cylinders are 
driven by two independent Brushless DC motors (Maxon, 
250W), at a speed up to lOHz. The torque T exerted 
on the inner cylinder is measured with a HBM rotating 
torquemeter (T20WN, 2N.m). A photograph of the setup 
is presented in figure [TJ The optical measurements have 
been done with an outer glass cylinder, of radius 121 ± 
0.25mm (r; = 0.909). 
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FIG. 1 Experimental setup, with the rotating torquemeter 
(upper part of picture), the calibration grid displacement de- 
vice (on top of the upper plate), one of the two cameras (left 
side) and the light sheet arrangement (right side) . The second 
camera is further to the right. 

For a newtonian fluid of kinematic viscosity v, we use 
the set of parameters defined by DubruUe et al. (2005): 
a mean Reynolds number Re = 2/(1 — rf) \(riroij-iod/ 1^) — 
{riUJid/v)\ based on the shear and on the gap d; and a 
"Rotation number" Ro which is zero in case of perfect 
counterrotation (r^w^ — —roOJo)- At lOHz in counter- 
rotation, the typical shear rate is around 1400s~^ and 
_Re ~ 1.4 X 10^ for pure Water. 

The physical properties of the fluids we use are given 
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TABLE I Physical properties of the fluids at 20°C; density p 
(kg.m"'^), dynamic viscosity p, (Pa.s) and kinematic viscosity 
V (m^.s-^). 



P M ^ 

Water 1000 1.0 x 10"^ 1.0 x lO"*' 

Nal solution 1500 2.0 x 10"^ 1.3 x 10-** 
Oil 800 3.0 X 10~-^ 3.8 x lO"*' 



counterrotation, a turbulent regime is easily reached (let 
us remind that the scaling expected in the Taylor- vortices 
regime is G oc Re^^^). 

Moreover, Esser and Grossmann (1996) proposed a for- 
mula for the first instability threshold. Here it gives 
RedrjjRo) — 338. At first glance, this value seems to 
be consistent with our results (see figure [^b) • In forth- 
coming experiments, we will use glycerol to lower the 
Reynolds number in order to check this prediction. 



in table HI The Oil is the Shell Macron EDMllO; it 
is a low-viscosity paraffinic hydrocarbons mix. Its re- 
fractive index measured with an Abbe refractometer at 
20°C is 1.4445. The interfacial tension between Oil and 
Water is cr = 0.045N.m~^. To have a refractive index 
matched dispersion (Budwig, 1994), in addition to pure 
Water, we also use a Sodium Iodide solution as acqueous 
phase (Narrow et al., 2000). The concentration of Nal 
is 510g.L~^, and traces of Na2S203 are added to avoid 
yellowish coloration. The viscosity ratio is then close 
to unity and the density ratio is close to 2. In the fol- 
lowing, unless specified, dimensionless quantities are de- 
duced from Ti as unit length, l/tOi as unit time, and prf 
as unit mass. 



IV. RESULTS 

A. Single-phase flow 

1. Torque scaling 

We first study the turbulent Taylor-Couette flow in ex- 
act counterrotation, and characterize its turbulent state. 
Though the Taylor-Couette flow has been widely stud- 
ied, few experimental results and theories are available 
for Ro — (DubruUe et al., 2005). We report torque 
measurements on pure fluids in figure. [H They have been 
made in a first version of the experiment, with a free sur- 
face and a space of 10mm between the cylinders bottom 
ends. The torque due to the bottom part has been re- 
moved by studying different filling levels. The range of 
Reynolds numbers we cover is 2.0 x lO'^ ^ < 1.2 x 10^. 
In figure[2^, we plot the dimensional torque vs. speed for 
the three pure fluids. The torques behave non-linearly 
with speed, as expected for turbulent regimes, and the 
Oil and Water are very close, the higher viscosity of Oil 
balancing its lower density. We plot in figure [^b the 
dimensionless torque G — T/{pv'^L) normalised by the 
laminar dimensionless torque Giam — 27777/(1 — rff'Re vs. 
Re. The experimental data fall on the same curve, which 
is consistent with the formula proposed by DubruUe and 
Hersant (2002) of the form G = aRe^ l{ln{bRe^))^'^/^\ 
Here, a non-linear fit gives a = 16.5 and b = 7 x 10"^. 
At the end of the range, this is close to a power-law be- 
haviour in G oc i?e^'^^. This scaling indicates that in 
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FIG. 2 (a): Torque iis. rotation frequency in exact counterro- 
tation for Oil (A), Water (o) and Nal (*). (b): Dimensionless 
torque G normalised by laminar torque Giam vs. Re. Solid 
line is the fit proposed by DubruUe and Hersant (2002). Rec 
is the expected value for the first instability threshold, ac- 
cording to the theory of Esser and Grossmann (1996). 



2. Stereoscopic PIV measurements 

We measure the three components of the velocity in 
one plane by stereoscopic PIV. The plane, illuminated 
by a double-pulsed Nd:YAg laser is normal to the mean 
flow, i.e. vertical (see figure [T]): the in-plane components 
are the radial (u) and axial (v) velocities, while the out- 
of-plane component is the azimuthal component (w). It is 
imaged using two double-frame PCO-cameras, mounted 
on Scheimpflug adaptors, located at each side of the light 
sheet, with an angle of 60° in the air. The flow is seeded 
with tracer particles (8 micrometer Sphericel). The field 
of view measures 11mm x 22mm, corresponding to a 
resolution of 300 x 1024 pixels. 

Special care have been taken concerning the calibra- 
tion procedure on which the third azimuthal compo- 
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nent rely. We use a very thin transparent sheet with 
crosses printed on it, with a weight at its bottom, and 
attached to a rotating and translating micro-traverse (see 
figure [T|). It is first put into the light sheet and traversed 
perpendicularly to it. Typically five pictures are taken 
with steps of 0.5mm. The raw PIV-images are processed 
with PivWare (Westerweel, 1993), with a last interroga- 
tion area of 32 x 32 pixels with 50% overlap, and nor- 
malised median filtering as post-processing. Then the 
vector mapping and the third component reconstruction 
are made with Matlab. The mapping function is a third- 
order polynomial, and interpolations are bilinear. 
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FIG. 3 Velocity profiles measured by stereoscopic PIV. Outer 
cylinder is at rest, Re ~ 90. Solid line: measured mean az- 
imuthal velocity. Dotted line: theoretical profile. Dashed 
line: fit of the form w = ar + h/r. The radial component u 
which should be zero is also shown. 

In order to check the quality of the technique, we 
measured the basic laminar flow with the outer cylin- 
der at rest, for {Ro = (1 - ?/) ; Re ~ 90}. The 
measurements are done in 86% glycerol-water mixture. 
We plot in figure [3] the velocity radial profiles, averaged 
over 100 recordings. The fiow is stationary and lami- 
nar. It is in good agreement with the theoretical profile: 
w{r) = {rj^/{'if — l))ijJi r + (1/(1 — rf'))ujirf 1/r, except 
near the inner and outer walls, where some refiections 
disturb the measurement. The error in the in-plane ve- 
locity is roughly constant, around 0.015. 

We then measured the counterrotating flow, at Re = 
1.4 X 10''. The results are presented in flgure [H The 
measurements are triggered on the outer cylinder posi- 
tion, and are averaged over 500 images. In the counter- 
rotating case, for this large gap ratio and at this value of 
the Reynolds number, the instantaneous velocity fleld is 
really desorganised and does not contain obvious struc- 
tures like Taylor- vortices, in contrast with other situa- 
tions (Wang et al., 2005). No peaks are present in the 
time spectra, and there is no axial-dependency of the 
time-averaged velocity field. We thus average in the ax- 
ial direction the different radial profiles. The azimuthal 
velocity w in the bulk is low, i.e. its magnitude is below 
0.1 between 112 < r < 120 that is 75% of the gap width. 
The two other components are zero within 0.002. 

The fluctuations of the velocity are not isotropic. They 




FIG. 4 Velocity measured by stereoscopic PIV for the coun- 
terrotating case, at Re ~ 1.4 x 10*. (a): Dimensionless 
mean velocity profiles u (dotted line), v (dashed line), and 
w (solid line). The azimuthal profile is not measured at the 
walls where its magnitude is 1. (b): Dimensionless turbu- 
lent shear-stress {u'w') (solid line) and dimensionless viscous 
shear-stress computed from the mean velocity profile. The 
dash-dotted line connects the two known values of the shear- 
stress at the walls, deduced from torque measurements. 



are of the order of 0.05 for the radial and axial compo- 
nents, and of the order of 0.1 for the azimuthal com- 
ponents. We can also compute some terms in the an- 
gular momentum balance equation (Marie and Daviaud, 
2004). According to the torque measurements described 
in the previous paragraph, the shear-stress at the inner 
wafl is Ti — T/{2TrLrf) ~ 1.49Pa, i.e. Ti/{prfomegaf) = 
3.1 X 10"'^ in dimensionless form. The dimensionless 
shear-stress at the outer wall is 2.8 x 10~^. We report 
these two values in figure [Hd, together with the viscous 
shear-stress pr{duj/dr) (dashed line) and the turbulent 
shear stress p{u'w') (solid line). The viscous shear-stress 
is negligible in the bulk, and is dominant at the walls, 
where it should ensure all the transport of angular mo- 
mentum. Though the boundary layer is very badly re- 
solved, the rough estimate of the viscous stress at the 
walls — corresponding to the tails of the dashed line — is 
consistent with the expected values. The measured tur- 
bulent shear-stress is roughly constant in the bulk, and 
has the good order of magnitude. The corre lation coeffi- 
cient between u and w, defined as {u'w') / ^ {u'u') {w'w') 
is 0.4. The turbulent shear-stress represents 85% of the 
expected shear-stress. We conclude that this fiow is dom- 
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inated by turbulent transport, which is also consistent 
with the high scaling of the torque with Re (Remind 
that G cx i?ei-^5). 



B. Two-phase flow 

1. Qualitative pictures of the dispersion process 

We now use two immiscible fluids in the Taylor- 
Couette facility. 




FIG. 5 LIF pictures of the dispersion in a counterrotating 
flow at 2Hz. Dark region corresponds to the Oil, and bright 
regions to the Sodium Iodide solution. The outer cylinder is 
on the right side of each picture. 
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FIG. 6 (a): Torque vs. rotation frequency for pure Oil (A) 
and 33% Water in Oil (*). Speed is increased by steps of 
0.2Hz every 60s, up to 4Hz and then decreased back to 0. 
The points clouds are parametric plots from experiments (b) 
and (c). (b): From a two-layer system at rest, the speed is set 
to 3.5Hz. (c): From this speed, after 20 minutes, the speed is 
decreased to 3Hz. 



To look at the phases repartition, we perform visuali- 
sations with a Laser Induced Fluorescence technique. We 
add a small amount of Rhodamine B in the Refractive- 
index matched Sodium Iodide solution. The Rhodamine 
is sensitive to the green light of the Nd:YAg Laser 
(532nm), and emits a yellow fluorescent light. We use 
Melles-Griot low-pass filters, with a cutoff wavelength of 
550nm to only record the fluorescent light. Examples of 
visualisations are given in flgure [51 

Starting from fluids at rest, the system is organised 
into two layers. When the cylinders are rotating, there 
is a competition between centrifugal forces and grav- 
ity. A Froude number measuring the ratio of the cen- 
trifugal forces to buoyancy can be estimated as Fr — 
~ 0.9/. We indeed observe that for / < IHz 
the fluids stay separated, but that the interface is tilted, 
the heavier fluid being pushed outwards. The inter- 
face is very weavy and non stationary. The ratio of 
turbulent shear stress to buoyancy can be estimated as 

Fr = i-^f'^ ^ 0.23/, with = 0.075 being the 
typical rms velocity measured by PIV. A Weber number 
comparing the effects of turbulent stresses to interfacial 

tension can be defined as We — ~ 0.6/^. At a 

speed of 1.5Hz, some Oil droplets start to be entrained 
into the acqueous phase. Very few drops of acqueous 
phase into Oil are present. Increasing further the speed. 



more and more droplets of Oil are entraped into Nal, 
and there start to be acqueous drops into the Oil, mainly 
visible as films, i.e. drops into drops, as can be seen in 
figure [5l Then, there seems to be three regions in the flu- 
ids, from top to bottom: an Oil-continuous region with 
Nal drops, a foam of possibly multiple droplets, and a 
Nai-continuous region with a lot of Oil droplets. Finally, 
for / > 3.5IIz, the fluids are fully dispersed. 

2. Evolution of the torque for dispersed phases 

One of the open question we adress is whether the tur- 
bulent dissipation is changed with respect to a monopha- 
sic case, and in that case, how much is that change. There 
are models for the effective viscosity of emulsions (Pal, 
2001), but what does happen in a turbulent dispersion 
where molecular viscosity is usually negligible ? 

We plot in flgure [6^ the torque vs. the frequency for 
a Water and Oil system. The volumic fraction of Water 
is 33%. The experiment have been automatised and we 
proceed as follows for this experiment: starting from rest 
we record the torque during 60s and then increase the 
speed by 0.2Hz, record one minute and so forth until 4Hz, 
which is reached after 20 minutes. We then decrease the 
speed until 0, by steps of 0.2Hz. 
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The torque for the Water and Oil system (★) is at low 
speed very close to the single-phase case (plotted in the 
same figure for pure Oil (A)). This is consistent with 
a sum of two layers weakly interacting (small viscosity 
difference), which basically give the same dimensionnal 
torque as discussed in figure [21 Increasing the speed, 
the torque shows up a fast increase around / = 3Hz 
(Fig. [5^) , corresponding to the point at which the liquids 
start to fully disperse. At the speed of 4Hz, the torque 
per unit mass of fluid to stirr is 0.27N.m.kg~^ whereas 
it is 0.15N.m.kg^^ for pure Oil. This corresponds to an 
increase of 70%. 

When decreasing the speed, there seems to be hystere- 
sis in the system. The torques on the return path are 
slightly higher than during the increase. The end of the 
hysteresis loop seems to be at / ~ 2Hz. Indeed, visual ob- 
servation confirms that the fiuids are still dispersed, even 
below 2Hz. Similar behaviours are present for other con- 
centrations we studied (see figure [7]), with less hysteresis 
for the lowest oil concentration. To better understand 
this hysteresis, we performed the following long-time ex- 
periments. Starting from a two-layer system at rest, we 
increase quickly the rotation up to / = 3.5Hz and record 
then the torque (figurellb). The torque gradualy increase 
during the 600s of the experiments. The final stationary 
state is not reached before the end of these 600s. We plot 
also this experiment in figure [6K as a parametric points 
cloud with time as parameter. One can notice that the 
torque at the beginning is the one for Oil, and gradually 
reaches the upper branch. 

The complementary experiment (figure [Ht) consists of 
lowering the speed to 2.8IIz after 20 minutes. It shows 
a slow decrease and the final state is not on the lower 
branch (figure [6^) : the two-layers system is not the one 
stable above 2.8Hz. The time scales involved in the dis- 
persion and separation processes are thus very long, the 
latter being even longer (figure [6)3-c) . Complicated be- 
haviours including periodic oscillations (with a period of 
20 min) have been observed in between and requires fur- 
ther investigation. 

We report in figure [7] the evolution with the volume 
fraction of the acqueous phase of the torque per unit mass 
of fluid to stirr at 4Hz. The data have been normalised 
by that of pure Oil (0.15N.m.kg~^). This quantity is 
maximum for a volume fraction of the acqueous phase 
around 33% and reaches a factor 2. Similar behaviours 
have been evidenced recently in liquid-liquid turbulent 
pipe flow experiments (Piela et al., 2006). 



V. CONCLUSION AND PERSPECTIVES 

We have built a facility with well-controlled parameters 
to study liquid-liquid dispersions in a turbulent flow. In 
a first step, the single-phase fiow has been characterized 
in terms of scaling of the dissipation, and also in terms 
of velocity and Reynolds stress profiles. We measured 
the turbulent fluctuations intensity who are a key feature 
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FIG. 7 Massic torque at 4Hz, normalised by the pure Oil case 
vs. volume fraction of Water (o) and Nal (★). 




FIG. 8 (a): Torque vs. rotation frequency for pure Water (o) 
and 50% Water in Oil (*). (b): Dimensionless Torque G vs. 
Re, molecular value for pure Water (o), and effective values 
proposed by Pal (2001) for 50% Water in Oil (★). 



for the dispersion mechanisms. Further investigations 
on the scaling and threshold for flrst instability for the 
counterrotating monophasic flow will be performed, on a 
wider Reynolds number range, using glycerol. 

A global measurement such as in-line torque measure- 
ment is able to characterize the state of the system in 
terms of whether it is dispersed or not. Data presented 
here for the torques in the dispersions were still limited by 
old motors with low power. Further work have to be done 
to explore the scaling of the torque for turbulent disper- 
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sions. We plot in figure [8] a first result obtained for 50% 
Water in Oil. The torque in the dispersed state seems 
to follow a nice scaling at higher velocities. In order to 
check the idea of a proper " effective molecular viscosity" , 
we try the one proposed by Pal (2001): for viscosity ratio 
3 (oil drops in water) and for fraction 0.5, the effective 
viscosity is 5.9 x 10~^m^.s~^. It seems then possible to 
collapse the points onto a single curve when one uses this 
effective viscosity to compute the Reynolds number and 
to normalise the torque in the dispersed state (figure[5jD). 
Notice that in fact we do not have enough points at low 
values of Reynolds numbers for newtonian fluids to really 
compare this rescaling with consistent data. 

In order to understand the interaction between the dif- 
ferent scales, we are now working on edge detection al- 
gorithm to study the drop size distribution with the LIF 
images. 
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